Spin and orbital effects of Cooper pairs coupled to a single magnetic impurity 
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The Kondo effect strongly depends on spin and orbital degrees of freedom of unconventional su- 
perconductivity. We focus on the Kondo effect in uniformly gapped superconducting systems (the 
two-dimensional p x + ip y -wave and d x 2_ y 2 + id xy -wa,ve superconductors here) to compare the mag- 
netic properties of the spin-triplet and spin-singlet Cooper pairs. The difference appears when both 
of the paired electrons couple to a local spin directly. For the p x + ip y -wave, the ground state is 
always a spin doublet for a Simp = 1/2 local spin, and it is always a spin singlet for Si mp = 1. The 
latter is due to uniaxial spin anisotropy of the triplet Cooper pair. For the d x 2_ y 2 + id xy -wa.ve, the 
interchange of ground states occurs, which resembles a competition between the Kondo effect and 
the superconducting energy gap in s-wave superconductors. Thus the internal degrees of freedom 
of Cooper pairs give a variety to the Kondo effect. 



I. INTRODUCTION 

The problem of magnetic impurities m .superconduc- 
tors has been investigated for a long timeJjtl This Kondo 
problem was solved at least for the case of a standard 
BCS (s-wave) superconductor. Since the superconduct- 
ing energy gap A competes with the Kondo tempera- 
ture Ik, the Cooper pairs are partly destroyed and a 
local spin tries to couple with the conduction electrons 
antiferromagnetically to form a local singlet (Kondo sin- 
glet). As a consequence, the Kondo singlet is stabilized 
for Tk/A 3> 1, while the local spin is almost free from 
the Cooper pairs for Tk/A <C 1. This shows that a 
singlet ground state can be realized for large Tk even 
if the energy gap exists. An interesting question arises 
in mind when we introduce spin and orbital degrees of 
freedom into Cooper pairs: is the Kondo singlet still sta- 
ble for large values of Tk/A in unconventional supercon- 
ductors? In general, the unconventional superconductiv- 
ity is characterized by a spin and momentum dependent 
order parameter A (7(7 /(fe). The k dependence is classi- 
fied into several angular momenta (s, p, d, /) or these 
combinations. A gapless superconductor expressed by 
An(k^—ky) is observed in the high-T c cuprates.El A possi- 
bility of spin-triplet superconductor is expected for some 
heavy fermion compounds and ruthernates: the simplest 
expression of the k dependence (orbital part) of the order 
parameter is given by (k x + iky), which is proposed for 
time-revepaLbreaking superconductors, UGe2cJ'til and 
Si'2Ru04.E3tj In the experimental side, much effort has 
been made to detect evidence of the Kondo effect below 
the superconducting transition temperature in high-T c 
cuprates. A recent nuclear magnetic resonance experi- 
ment reported that a local magnetic susceptibility shows 
reduction of Kondo screening due to opening of the su- 
perconducting energy gapo 

In this paper, we focus on two simple unconventional 
superconducting states which have a full energy gap and 
whose angular momenta are good quantum numbers, as- 
suming that a local spin is located at the center of two- 



dimensional coordinates. One is a p x + ip y -wave state 
whose order parameter is given by A^(k x + iky), and 
the other is a d x 2_ y 2 + id xy -wave one represented by 
Af|([fc x 2 — k y 2] + ik x ky). Since the density of quasipar- 
ticle states for both cases has the same energy gap with 
the s-wave, we can extract an orbital effect of Cooper 
pairs which is never seen in the s-wave superconductiv- 
ity. We can clarify a spin effect in the Kondo effect as 
well to compare the spin-triplet pair with the spin-singlet 
one. We investigate a new type of Kondo effect due to 
a Si m p = 1/2 or S lmp = 1 impurity (Si mp is a size of 
the local spin) in these superconducting states, us ing the 
numerical renormalization group (NRG) methodML3 In 
our NRG model, it is found that the quasiparticles of the 
superconducting states have two channels which are cou- 
pled to each other only at the impurity. The NRG tech- 
nique is very powerful to study such complicated Kondo 
problems. 

In our previous study on the_Kondo effect in the uncon- 
ventional superconductors J13c3 a short-range scattering 
object was considered for the impurity. The atomic or- 
bital of the impurity is restricted to the s-orbital denoted 
by the I — angular momentum. The angular momen- 
tum of the Cooper pair is I = 1 for the p x + ip y -wave 
and I — 2 for the d x 2_ y 2 + id xy -wa,ve. In metallic states, 
I = is the only relevant orbital in the Kondo effect. For 
the p x + ipy-wave, the Cooper pair consists of I = and 
I = 1 electrons and involves the I = 1 electrons in the 
Kondo effect. In the same manner, I = 2 electrons be- 
come relevant for the d x 2_ y 2 + id xy -wewe. This leads to a 
doublet ground state and the situation is completely dif- 
ferent from the s-wave superconducting case. However, 
in this case, we did not find any difference between the 
Px + ipy-wave and d x 2_ y 2 + id xy -w&ve cases. 

In the present study, we take into account the scatter- 
ings of the I = 1 (I = 2) electrons on the impurity site as 
well as I = for the p x + ip y -w&ve (d x 2_ y 2 + id xy -wave). 
Introducing the two exchange couplings, we can distin- 
guish the spin-triplet Cooper pair from the spin-singlet 
one in the Kondo effect. The difference appears in the 
form of the effective exchange interaction in our Kondo 
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model. The p x + ip y -wave generates spin anisotropy, and 
the ground state is always a spin doublet for Si mp = 1/2. 



For Si, 



1, the triplet state of the local spin is split into 



a singlet and a doublet, and the ground state is always a 
spin singlet. For the d x 2_ y 2 -Hcfcy-wave, the ground state 
changes from a spin doublet to a particle-hole doublet as 
Tk increases for Simp = 1/2. Similarly, the interchange 
of triplet and singlet ground states occurs for Simp = 1. 

This paper is organized as follows. The NRG formu- 
lation of our Kondo problem is presented in Sec. II. In 
Sec. Ill, the difference between the two types of uncon- 
ventional superconductivity is discussed, based on the 
NRG results. The paper is closed with concluding re- 
marks in Sec. IV. 



II. MODEL 

First we summarize the derivation of the NRG Hamil- 
tonian for the Kondo effect in. the p x +ip y -w&ve (d x 2_ y 2 + 
id xy -w&ve) superconductor £j One can find that this type 
of superconducting order parameter is uniformly gapped 
since A.(k x + iky) for the p x + «p y -wave is expressed with 
Ae^fc , where is the angle of a wave vector k measured 
from the k x axis and A is a real function of |fe|. We note 
that this A is not a function of k itself in this case. 

Let us begin with the following Hamiltonian: 
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Here -f/kin, Ha, and H? mp represent the kinetic energy 
of conduction electrons, the BCS interaction, and the ef- 
fective exchange interaction at the impurity, respectively. 
The operator aj^ (afc CT ) represents creation (annihilation) 
of a conduction electron with a wave vector k and spin 
a. In H? mp , a is the Pauli matrix for the conduction 
electrons. The impurity is expressed by the local spin 
operator S. The exchange coupling J(< 0) is antiferro- 
magnetic. Since the orbital angular momentum is a good 
quantum number for the p x + ip y -wave Cooper pairs, we 
can simplify the Hamiltonian by expanding a^a with re- 
spect to the two-dimensional polar coordinate bases: 



(2.5) 



where Ji is the Z-th Bessel function, R is the radius of the 
two-dimensional system, and I is the ^-component of the 
orbital angular m ome ntum of the conduction electron. 
The Hamiltonian (pj) is then expressed as 
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In Eq. (2.8), we have assumed a short-range scattering 
impurity. It is sufficient to choose the I = and / = 1 
orbitals and neglect others for the p x + ip y -wetve case, 
since they are the only orbitals of conduction electrons 
that participate in the Kondo effect. The I = 1 electrons 
are coupled to the local spin via the superconducting or- 
der parameter. The most appropriate forms of the three 
parts of Hamiltonian ( |2.1[ ) for the NRG calculation are 
finally given by 
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A similar Hamiltonian for the d x 2_ y 2 + id xy -w&ve can be 
obtained if iA and I = 1 are replaced by — oA and / = 2, 
respectively, since the order parameter A([fe x 2 — k y z] + 
ik x k y ) is expressed with Ae l2( ^k. 

Next we present the derivation of the NRG Hamilto- 
nian from the Kondo model. For simplicity, we use the 
same values for both the superconducting cutoff and the 
band width, since it is assured that this does not alter 
the results.Q After applying the Wilson's logarithmic dis- 
cretization procedure to the conduction band,tj we ob- 
tain the following hopping type of Hamiltonian: 
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for the BCS pairing interaction in the p x + ip y -wa,ve su- 
perconductor. To diagonalize Ha, we use the following 
Bogoliubov transformation: 



and we obtain 
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By applying the following particle-hole transformation^ 
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we obtain the standard form of NRG Hamiltonian in a 
recursion relation 
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where p is the density of states at the Fermi energy and A 
is normalized by the conduction band width here. The 
above procedure can be applied to the d x 2_ y 2 + id xy - 
wave case, and the same form of the NRG Hamiltonian 
is obtained. 

In the above argument, we have restricted an orbital 
type of the scattering impurity to the s-wave (I = 0). 
In this case, both p x + ip y -w&ve and d x 2_ y 2 + id xy -w&ve 
superconducting stat es ex hibit t he sa me Kondo behavior 
resulting from Eqs. (2.18) and (2.19). The simplest way 
to distinguish the two types of pairings is that we add 
p-wave scattering (1 = 1) for the p x + ip y -w&ve case and 
d-wave scattering (I = 2) for the d x 2_ y 2 + id xy -w&ve. 
First, we discuss an extended Kondo Hamiltonian for the 
former, which is given by 

-frfmp = X/ ^ ' ^ct'C-^o/oOct/oOct _ JlfoiafoiT'), (2.22) 
era' 

where Jo and Ji are coupling constants for the I = and 
I = 1 electrons, respectively, and the fermion operator 
folo- is defined by 
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From Eqs. fl2.14D , ( |2.15[ ), and ( p7l7| ), we obtain 
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and transform to the NRG Hamiltonian for the 
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where S± = S x ± iS y , and the orbitally dependent ex- 
change couplings are given by 



J* T , = Jo + sgn(rr')Ji, 
J TT , = Jo - sgn(rr')Ji- 



(2.27) 



We can find that uniaxial anisotropy is introduced to the 
spin of the total system through the couplings: J* T , ^ 
J^ T , for finite values of J and J\. On the other hand, for 
the d r 2_„ 2 + d xy -wsve, the transformation corresponding 
to Eq. ( 2.25| ) is given by 
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for the I = 2 electrons. The Hq for the Z = and I = 2 



electrons is given by the same form with Eq. (2.26) except 
for the exchange couplings 



J* T , = J + sgn(rr')J2, 
Jtt' = Jo + sgn(rr')J2, 
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where J2 is a coupling constant between the impurity 
and I = 2 electrons. In this case, the total spin keeps its 
isotropy: J* T , = J^ T , is satisfied for all values of Jo and 
J 2- Thus, by introducing the I 7^ scattering into the ex- 
change interaction, we can see the difference between the 
triplet and singlet pairings which comes from the phases 
of the wave functions of the I ^ electrons. 



III. RESULTS 

A. Short-range scattering impurity 

First we summarize main results for a case of short- 
range scattering impurity with S- lmp = 1/2, the-details of 
which were discussed in our previous papers .EZHia Both 
Jo and A are relevant in the Kondo effect, growing as 
k N f 2 with the renormalization step N of the Hamilto- 
nian (2.18). The competition between the Kondo effect 
and the superconductivity is found in Xk/A dependence 
of the bound state energy. The ratio of this energy and 
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FIG. 1. Ji (< 0) dependence of the bound state energy 
levels for Simp = 1/2 in the p x + ip y -wave superconducting 
state, where J = -0.2 and A = 10" 5 (A = 1.5A). The 
energies are measured from those of the spin doublet ground 
state (S z = ±1/2). The two bound states are particle- hole 
doublets (Sz = 0). The circles and stars represent data for 
one bound state dominated by Jo and the other dominated 
by Ji, respectively. In the NRG calculation, the lowest-lying 
~500 states are kept at each renormalization step. 



the renormalized A approaches a constant value when A 
becomes large enough. In a logarithmic scale of Tk/A, 
the convergent value (-Eb/A)* measured from a spin- 
doublet ground state changes monotonically from a unity 
for Tk/A -C 1 to zero for T K /A > 1. This means that 
the ground state is always a spin doublet for the p x + ip y - 
wave superconductivity. However, the Kondo effect does 
occur since Tk is found to be an important energy scale. 
In fact, the local spin is shrunk as Tk/A increases. It is 
shown by the decrease of an effective ^-factor which is es- 
timated frora-Curie-law behavior of the local spin suscep- 
tibility xioc-li3 In the strong coupling limit (Tk/A 3> 1), 
the local spin is almost quenched by the I = electrons, 
while the I ^ electrons couple weakly with the local 
spin via the superconducting order parameter to gain the 
superconducting condensation energy. Therefore, a sin- 
glet ground state is not realized even if Tk is very large, 
which is completely different from the s-wave supercon- 
ducting case. For the latter, the interchange of singlet 
and doublet ground states occurs. 

In the above result, we did not see any difference be- 
tween the triplet (p x + ip y -w&ve) and singlet (d x 2_ y i + 
id xy -w&ve) pairing states. For both pairings, the orbital 
effect of the Cooper pairs is important in the Kondo ef- 
fect, while their spin does not affect the Kondo effect. 



B. Multiorbital scattering impurity 

Next we discuss a case of multiorbital scattering im- 
purity to distinguish the two types of superconductivity 
in the Kondo effect. Fo r thi s purpose, we have derived 
the NRG Hamiltonian ( |2.26| ). For A = and J = Ji 
(J2), tiifi model corresponds to the two-channel Kondo 
model .cj In this metallic case, the Kondo temperature 




FIG. 2. Tk/A dependence of the energies of the low- 
est-lying particle-hole doublet (S = 0) and the second low- 
est-lying spin doublet (S = 1/2) for Simp = 1/2 in the 
d x 2_ y 2 + id X y-wave superconducting state. The energies 
are measured from those of the lowest-lying spin doublet 
(S = 1/2). The circles (Jo = J2 = —0.2) and squares 
(Jo = J2 = —0.3) represent the data for the S = state. The 
plus (Jo — J2 ~ —0.2) and cross symbols (Jo = J2 = —0.3) 
represent the data for the second lowest-lying S — 1/2 state. 
The interchange of the S = 1/2 and S = ground states oc- 
curs at Tk/A ~ 0.3. In the NRG calculation, the lowest-lying 
~500 states are kept at each renormalization step. 



is given by Tk = | Jop\ exp(— 1/\ Jop\)- In the logarith- 
mic discretization procedure for NRG, JqP is effectively 
reduced by the factor g(A):Q g(3) = 0.910 is taken here. 



1. S in 



1/2 impurity case 



Let us begin with the Si mp = 1/2 impurity. For the 
Px + ipy-wave, the ground and the first excited states 
are always spin-doublet (S — 1/2) and particle-hole dou- 
blet (S = 0), respectively, which is same as the J\ = 
case discussed above. The Tk/A dependence of the first 
excited state energy is also similar, but it does not ap- 
proach zero for large values of Tk/A. Figure [l] shows 
the Ji dependence of the first and second excitation en- 
ergies when Jo and A are fixed at 0.2 and 10 -5 , respec- 
tively, which give T K /A = 13.8 for J x = 0. The first 
excited energy (Eb.i/A)* given in Fig. [I] increases grad- 
ually with increasing | Ji| for | Ji| < 0.2. This means that 
the effect of p-wave scattering does not destabilize the 
Ji = state immediately The second excitation energy 
level (-Eb,2/A)* is located above 1.0. When Ji is close to 
<T>, (Tb,i/A)* upturns suddenly and (£3,2/ A)* appears 
below 1.0 as the second bound state. Both energy lev- 
els cross just when Ji is equal to Jo- When J and J\ 
are exchanged, the same result is obtained. For A = 0, 
we obtain the non-Fermi liquid fixed point derived from 
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FIG. 3. J2 (< 0) dependence of the energy of the low- 
est-lying particle-hole doublet (S = 0) for Si mp = 1/2 in the 
d x 2_ y 2 + id X y-wa,ve supercinducting state, where Jo is fixed at 
—0.2 (A = 1.5A). The energies are measured from those of 
the spin doublet (S = 1/2). In the NRG calculation, the low- 
est-lying ~500 states are kept at each renormalization step. 

the two-channel Kondo mode0 for Jo = J±. Once the 
two couplings are anisotropic in this case, the non-Fermi 
liquid fixed point becomes unstable immediately, and a 
Fermi-liquid fixed point is stabilized by a single relevant 
channel. When A is finite for J± — Jo, the two indepen- 
dent channels are combined by A. The continuous change 
of the curve (£"b.i/A)* in Fig. [I] means that there is no 
discontinuous change of the low-energy spectrum in all 
the parameter space of Jo and Ji . The abrupt change of 
(Tb,i/A)* around Ji = Jo corresponds to the crossover 
from the Jo dominant spectrum to the Ji dominant one. 

Contrary to the above p x + ip y -w&ve case, the particle- 
hole doublet (S = 0) can be the ground state for the 
d x 2_ y 2 + id xy -w&ve when Tk is large enough. Figure ^ 
shows the Ik/ A dependence of the particle- hole doublet 
(S = 0) energy (E/A)* measured from the spin doublet 
(S = 1/2) for Jo = J 2- The interchange of the ground 
states is very similar to a S- lmp = 1/2 impurity in the s- 
wave superconductor where only I = orbital is involved. 
We note that the particle-hole doublet (S = 0) ground 
state is not smoothly connected to the A ^ILlimit where 
a spin-doublet ground state is realized ESE3 As IJ2I is 
decreased from | J2 1 = | Jo | , the spin-doublet state is more 
stable as shown in Fig. || In other words, the coupling J 2 
lowers the energy of the Kondo singlet formed against the 
superconducting condensation energy. Thus the role of 
the coupling J2 for the d x 2_ y 2 + id xy -wave is completely 
different from that of the coupling Ji for the p x + ip y - 
wave when their magnitude is close to | Jo | . In Fig. [jj the 
energy of the first excited state (S = 0) increases as | Ji| 
(< I Jo |) becomes larger, implying that the Kondo singlet 
becomes unstable. We conclude that the roles of I ^ 




FIG. 4. Xk/A dependence of the bound state energy for 
Simp = 1 in the p x + ip y -wa,ve supercinducting state. The en- 
ergies are measured from those of the spin singlet ground state 
(S z — 0). The bound state is a spin doublet (S z = ±1). The 
circles and squares represent the data for Jo = Ji = —0.15 
and for Jo = Ji = —0.2, respectively. In the NRG calculation, 
the lowest-lying ~800 states are kept at each renormalization 
step. 



exchange couplings for the spin-singlet and spin-triplet 
superconductors are opposite to each other for Si mp = 
1/2. 



2. S in 



1 impurity case 



When T K /A < 1 is satisfied for S [mp = 1/2, the 
spin-doublet is the ground state in both superconduct- 
ing states. We can see a more drastic difference between 
them for the Si mp = 1 case discussed below. The ex- 
change couplings of I = and / = 1 electrons with the 
local spin breaks the spherical symmetry of the total spin 
in the p x + ip y -wsve superconducting state, and they lift 
the degeneracy of spin triplets. The anisotropy is uni- 
axial and a spin singlet is stabilized against a doublet 
in this case. As a consequence, a singlet ground state 
is realized even for Tk/A -C 1. Figure [| shows the 
Tk/A dependence of the energy of the first excited dou- 
blet state (S z = ±1) measured from the singlet ground 
state (S z = 0) for J 1 = J . This singlet for T K /A < 1 is 
not related to a Kondo singlet but to the spin anisotropy 
introduced by the exchange coupling J\ in addition to Jo. 
When Tk/A is large enough, there is no impurity bound 
state below the energy gap. In this sense, the local spin 
is compensated by the Kondo effect. The ground state is 
always a spin singlet except for either Jo = or J\ = 
case where it is a spin triplet over all the values of Tk/A. 

On the other hand, for the S\ mp — 1 impurity in 
the d x 2_ y 2 + id xy -w&ve case, we find the interchange of 
spin-triplet and spin-singlet ground states as found for 
Simp = 1/2. An energy level of a quartet state (degen- 
erate S = 1/2 doublet states) appears between those of 
the singlet and triplet states as shown in Fig. ||. The 
quartet (S = 1/2) cannot be the ground state for all the 
values of Tk/A. As | J2I is decreased from \Jq\, the spin- 
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FIG. 5. Tk/A dependence of the singlet-state (5 — 
state) and quartet-state (twofold degenerate S = 1/2 state) 



FIG. 6. Ji (< 0) dependence of the singlet-state (5 = 
state) and quartet-state (twofold degenerate 5 = 1/2 state) 



energies for S- w 



1 in the d x 2_ 2 + id xy -wsve supercinduct- energies for 5i mp = 1 in the d x 2_ y 2 + id xy -wa,ve supercinduct- 



ing state. The energies are measured from those of the spin 
triplet (5 = 1). The circles (Jo = J2 = —0.15) and squares 
(Jo = J2 = —0.2) represent the energies of the singlet (5 = 0). 
The plus (Jo = J 2 = -0.15) and cross (J = J 2 = -0.2) 
symbols represent those of the quartet (5 = 1/2). The in- 
terchange of the 5 = 1 and 5 = ground states occurs at 
T K /A ~ 0.05. In the NRG calculation, the lowest-lying ~800 
states are kept at each renormalization step. 



triplet state is stabilized as the spin-doublet is stabilized 
for Simp = 1/2. This is shown in Fig. |[ Contrary to the 
Simp = 1/2 case, the difference between the two super- 
conducting states is remarkable for Tk/A -C 1. This is 
due to the spin anisotropy introduced by the two-orbital 
couplings of the Kondo effect in the p x + ip y -wave su- 
perconducting state, while the spin of the total system is 



isotropic in the d x 2_ x 



ing state. The fixed parameters are Jo = —0.15 and A = 10~ J 
(A = 1.5A). The energies are measured from those of the spin 
triplet (5 = 1). The circles and stars represent the data for 
the singlet (5 = 0) and quartet (5 = 1/2) states, respectively. 
In the NRG calculation, the lowest-lying ~800 states are kept 
at each renormalization step. 



pie with the local spin. On the other hand, Xi oc (T — ► 0) 
shows a constant value and Ximp(T — * 0) vanishes when 
a spin-singlet ground state is realized for S imp = 1 in the 
Px + ipy-w&ve superconducting state and for small values 
of A in the d x 2_ y 2 + id xy -wa,ve. For the p x + ip y -wave, 
Xioc decreases with the increase of A, implying the Van 
Vleck like behavior. 



IV. CONCLUSION 



C. Impurity susceptibility 

Finally we mention the different low-temperature be- 
havior between two types of magnetic susceptibility, Ximp 
and xioc in the above cases. The former is the impu- 
rity susceptibility defined by Ximp — ix ~ Xo) where 
X (xo) i s the magnetic susceptibility of the conduction 
electron system including (in the absence of) the local 
spin. The latter is the local spin susceptibility defined by 
limh^o(M z /h) where M z — <?/ib < S z > is the magneti- 
zation of the local spin and h is a magnetic field coupled 
to the local spin. For S- lmp = 1/2 in the p x + ip y -wave, 
the ground state is a spin doublet for all the values of 
Tk/A, as we have discussed. In this case, both Xi mp 
and Xioc shows the Curie law. In zero temperature limit 
(T — ► 0), Txioc decreases with the increase of Tk/A, 
while Tximp always reaches a constant S , i mp (S , i mp + l)/3. 
Such difference indicates that the local spin is shrunk as 
the Kondo effect is stronger and the Cooper pairs cou- 



We have studied a S- lmp — 1/2 or S lmp = 1 impurity 
coupled to the p x + ip y -wa,ve or d x 2_ y 2 +id xy -wave super- 
conducting state in order to find the different behavior 
between the spin-triplet and spin-singlet Cooper pairs in 
the Kondo effect. The Kondo effect in both types of su- 
perconductors can be described by the same model for a 
short-range scattering impurity with only I = electrons. 
In this case, only one of the paired electrons couples to 
the local spin directly. It means that this type of effec- 
tive exchange interaction does not depend on the spin of 
the Cooper pair. To see spin effects of Cooper pairs, it is 
necessary to add the I ^ orbital of the paired electrons 
to the effective exchange interaction. We consider I = 1 
for the p x +ip y -w&ve and / = 2 for the d x 2_ y 2 -Hd^-wave. 
The two-orbital exchange interaction generates uniaxial 
spin anisotropy for the p x + ip y -wave since the spin of 
the Cooper pair is equal to one and is aligned to the x-y 
basal plain. On the other hand, spin isotropy is conserved 
for the d x 2_ y 2 + id xy -wave, since the singlet Cooper pair 



G 



has no spin. The two types of Cooper pairs cause the 
different Kondo behavior as follows: 



(1) 



1/2 local spin 



The ground state is a spin doublet in all the parame- 
ter region of Jo(< 0), Ji(< 0), and A for the p x + ip y - 
wave. For the d x 2_ y 2 + id xy -~wave, the interchange of the 
spin-doublet (for large A) and spin-singlet (for small A) 
ground states is found. 



(2) S h 



1 local spin 



The ground state is a singlet represented by S z = 
in all the parameter region as mentioned above for the 
Px + ipy-wave. It is split from a spin triplet due to the 
uniaxial spin anisotropy. For the d x 2_ y 2 + id xy -w&ve, the 
ground state interchanges between the spin-triplet (for 
large A) and spin-singlet (for small A) states. 

We have found that physical quantities such as bound 
state energies are scaled by Tk/A, which is a common 
feature for such competition between the Kondo effect 
and the energy gap. 

In order to apply our Kondo model to real systems, 
we must start from an extended Anderson model with 
orbital degeneracy. If we identify the orbital denoted by I 
with one of the atomic orbitals of the impurity, we have to 
take I = —1 (I = —2) into account as well as I = 1 (I = 2) 
for the p x + ip y -~w&ve (d x ?_ y i + id xy -wave). Since the 
p x +ip y -w&ve superconducting order parameter makes I = 
— 1 couple with 1 = 2, we need to include the four orbitals 
(I = —1,0,1,2) of conduction electrons in the present 
model. In the same manner, four orbitals (/ = —2, 0, 2, 4) 
of electrons are necessary for the d x 2_ y 2 +id xy -wave case. 
As for another point, the Simp = 1 impurity is actually 
obtained for the strong Hund coupling limit. If the Hund 
coupling is weak, effective orbital exchange interactions 
should be present. Since the conduction-electron part 
including the BCS interaction term is tran sform ed to the 
hopping type of NRG Hamiltonian in Eq. (2.18) first, the 
orbital exchange interactions break the conservation of 
the total number of quasiparticles and the z component 
of the total spin. This makes the mumerical calculation 
much harder. Thus the model becomes more complicated 
for the orbitally degenerate impurity. Neverthless, the 
present form of our model can describe a local moment 
in superconducting states if each orbital can be regarded 
as an independent branch of a conduction band instead of 
an atomic orbital. It is interesting to search for a realistic 
case where such extention of our study is possible. 
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